Textbook Examples of Recursion 

by Donald E. Knuth 

Abstract. We discuss properties of recursive schemas related to McCarthy's "91 function" 
and to Takeuchi's triple recursion. Several theorems are proposed as interesting candidates 
for machine verification, and some intriguing open questions are raised. 
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John McCarthy and Ikuo Takeuchi introduced interesting recurrence equations as they were 
exploring the properties of recursive programs. McCarthy's function [7] 



f(x) = if x > 100 then x - 10 else f(f(x + 11)) 

wo: 

has become known as the "91 function," since it turns out that f(x) = 91 for all x < 101. Takeuchi's 
function [13] is a triple recursion 



t(x, y, z) = if x < y then y else t(t(x — 1, y, z), t{y — 1, z, x), t(z — 1, x, y)) , 

u 

which has proved useful for benchmark testing of Lisp systems because the recursion terminates only 
after the definition has been expanded a large number of times (assuming that previously computed 
values are not remembered). Neither of these functions is of practical importance, because no 
reasonable programmer would ever want to carry out such recursive computations on a realistic 
problem. Yet both functions are quite instructive because they illustrate important problems and 
techniques that arise when we consider the task of verifying computer programs formally. Therefore 
they make excellent examples for textbooks that discuss recursion. 

The purpose of this paper is to obtain new information about f(x) and t(x,y,z) and about 
, several closely related functions. Several of the theorems proved below should provide good test 

material for automated verification systems. A few open problems are stated, illustrating the fact 
that extremely simple recursions can lead to quite difficult questions. 

1. The 91 function. It is appropriate to begin by studying the 91 function, because 1991 is 
the year of John McCarthy's 64th birthday (and because a computer scientist's most significant 
birthday is the 64th). McCarthy originally wrote down the definition of f(x), as shown above, 
because he wanted to study a simple recursion whose properties could not be deduced by ordinary 
mathematical induction. After studying the definition, he was pleasantly surprised to discover that 
it had the totally unexpected "91 property." 

The 91 function certainly belongs to the set of significant textbook examples, because it is 
mentioned on at least 14 pages of Zohar Manna's well known text, Mathematical Theory of Com- 
putation [6]. The first published discussions of the function appeared in 1970 [7, 9], after it had 
been investigated extensively at Stanford's AI laboratory during 1968 [8]. 

Instead of using McCarthy's original definition, let's change the specifications a bit and consider 
the following function (see [6, Problem 5-8]): 

fix) = if x > 100 then x - 10 else f 91 (x + 901) 
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where f 91 (y) stands for /(/(••• (f{y)) ■••))' * ne 91-times-repeated application of /. According 
to this new definition, we have 

/(91) = / 91 (992) = / 90 (982) = • • • = / 2 (102) = /(92) . 

And a similar derivation shows that if 90 < x < 100 we have 

f(x) = f\x + 901) = • • • = f\x + 11) = f(x + 1) ; 

hence 

/(90) = /(91) = ■ ■ ■ = /(100) = /(101) . 
And /(101) = 91, so we have proved in particular that 

/(91) = 91. 

Now let's evaluate f(x) when x is extremely small, say x = — 10 6 . We have 

/(-1000000) = / 91 (-999099) = / 181 (-998198) = • • • = / 99811 (-791) 

= / 99901 (110) 
= / 99900 (100) 

and we know that /(100) = 91; hence 

/9 9900 (100) = ^99899(9^ = /99898 (91) = . . . = /(Q1) = gl _ 

In general, if x is any integer < 100, let m be the smallest integer such that x + 901m > 100, and 
let n be the smallest integer such that x + 901m — lOn < 100. Then m > 1, n < 91, and 

f( X ) = f 1+90m ( X + 901m) = f+^m-n^ + 9Qlm _ 1Qn ) = j90m-«( 9] ) ^ 

where the last step follows since 91 < x + 901m — lOn < 100. We conclude that f(x) = 91. (The 
final step is omitted if m = 1 and n = 91; that case occurs if and only if x = 100.) 

How many iterations are needed to compute f{x) by this definition, if we continue to apply 
the recurrence even when evaluating f(x) for values of x that have already been considered? Let 
F(x) count the number of times that the test 'if x > 100' is performed; then we have 

F(x) = if x > 100 then 1 else 1 + F(x + 901) + F(f(x + 901)) + 

+ F(f 2 (x + 901)) + • • • + F(f 90 (x + 901)) . 

(This is a special case of the general notion of a derived function, which is always jointly recursive 
with the function from which it has been derived; see McCarthy and Talcott [11].) A bit of 
experimentation reveals that F(x) also reduces to a simple function: 
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Lemma 1. F(x) = if x > 100 then 1 else 9192 - 91s. 
Proof. If x < 100 we have 

90 

F(x) - F(x + 1) = ^2(F(f k (x + 901)) - F(f k {x + 902))) . 

fc=0 

Now if x + 901 < 100, the sum reduces to F(x + 901) — F(x + 902), because the terms for k > 
are F(f k (x + 901)) - F(f k (x + 902)) = F(91) - F(91) = 0. In this case we let x' = x + 901. 
On the other hand if x + 901 > 100, let n be minimal such that x + 901 — lOn < 100. Then 
1 < n < 90, and F(f k (x + 901)) - F(f k (x + 902)) = F(91) - F(91) = for all k > n. We also 
have F(f n (x + 901)) - F(f n (x + 902)) = F(x + 901 - lOn) - F(x + 902 - lOn); and F(f k (x + 
902)) - F(f k (x + 901)) = 1 - 1 = for all jfe < n. In this case we let x' = x + 901 - lOn. In both 
cases we have found an x' such that 

F(x + 1) - = + 1) - F(x') , x < x' < 100 . 

The proof is therefore complete by induction on 101 — x if we simply verify that F(100) — F(101) = 
91. □ 

The 91 function suggests that we consider the more general recursive scheme 

f(x) = if x > a then x — b else f c (x + d) , 
where a is an arbitrary real number, b and d are positive reals, and c is a positive integer. 

Theorem 1. The generalized 91 recursion with parameters (a,b,c,d) defines a total function on 
the integers if and only if (c—\)b<d. In such a case the values of f(x) also obey the much simpler 
recurrence 

f(x) = if x > a then x — b else f(x + d — (c — 1) b) . 

Proof. It is not difficult to show that any function satisfying the generalized 91 recursion for c > 1 
must also satisfy 

f(x) = if x > a then x — b else f c ~ l (x + d — b) . 
For if x < a, let n be minimal such that x + nd > a. Then 

f°(x + d) = /c+(«-l)(c-l) ( x + nd ) = fn(c-l) ( x + nd _fy. 
f c -\x + d-b) = fc- l +(n-l)(c-l) ( x + d _ b + ( n _ y ^ . 

hence f c (x + d) = f c ~ 1 (x + d— b), as desired. 

To complete the proof, we use induction on c; and we also need to characterize the parameter 
settings that cause the given recursive definition to terminate for all x. 

If (c — 1) b > d, the expansion of f(x) will not terminate when a — b < x < a. For if n is 
minimum such that x + d — nb < a, we have n < c — 1, and 

f{x) = f c (x + d) = ■ ■ ■ = f c - n (x + d-nb). 
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Now c — n > 0, and a — b < x + d — nb < a, so this will go on and on. 

On the other hand, we can show that no looping will occur if (c — 1) b < d. Suppose first that 
x > a — b. If x > a, obviously f(x) = x — b. Otherwise we have x + d > x + (c — 1) b > a + (c — 2) b, 
hence 

f(x) = f c (x + d) = ■ ■ ■ = f 2 (x + d - (c - 2) b) = f(x + d - (c - 1) 6) . 
Let A = d — (c — 1) 6, and let m be minimum such that x + m A > a; then 

/(x) = /(x + A) = ■ • • = f(x + m A) = x + m A - b . 

Thus, the expansion of f(x) terminates with a value > a — b whenever x > a — b. 

Finally, if x < a — b and if m is minimal such that x + md > a — b, the expansion of 

f{x) = f l+m( > c - l Xx + md) 

terminates, because we can peel off the /'s one by one. □ 

When the generalized 91 function is total, we can express it in "closed form" as 

f(x) = if x > a then x — b 

else a + d — cb — ((a — x) mod (d — (c — 1) 6)) . 

The special case c = 2 of Theorem 1 was first proved by Manna and Pnueli [8] . 
Open Problem 1. Prove Theorem 1 by computer. □ 

2. The Takeuchi function. Now we turn to the more complex recurrence 

t(x, y, z) = if x < y then y else t(t(x — 1, y, z), t(y — 1, z, x), t(z — 1, x, y)) . 

John McCarthy observed in unpublished notes [10] that this function can be described more simply 
as 

t(x, y, z) = if x < y then y else if y < z then z else x . 

If we assume termination, the latter function satisfies Takeuchi's recurrence, so it must be identical 
with the former function. 

John had just returned from a conference in Kyoto, and his notes [10] began with a brief 
comment about the history of this function and its motivation: 

Ikuo Takeuchi (1978) of the Electrical Communication Laboratory of Nippon Telephone 
and Telegraph Co. (Japan's Bell Labs) devised a recursive function program for comparing 
the speeds of LISP systems. It can be made to run a long time without generating large 
numbers or using much stack. 

(Incidentally, I believe [10] was John's first experiment with the use of T^X, a computer typesetting 
system that I was developing while sitting in the office next to his. Without his generous provision 
of computing and printing facilities, T^X would never have existed.) 
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At about the same time, John coerced the FOL proof-checking system to construct a 50-step 
proof that A(x,y,z) has the simple form stated above [10]. This experiment suggested several 
improvements to FOL. 

If we fully expand the definition of t(x, y, z) whenever x > y, the proof of termination seems 
to be nontrivial, because there is no obvious way to impose an order on the set of all arguments 
(x, y, z) in such a way that no infinitely long dependency chains exist. We shall prove termination 
as a byproduct of a more general investigation of the total running time needed to evaluate t(x, y, z) 
by repeated application of the definition. 

Let T(x, y, z) be the number of times the else clause is invoked when t(x, y, z) is evaluated by 
ordinary Lisp recursion. Then 

T(x, y, z) = if x < y then 

else 1 + T(x- 1, y, z) + T(y - 1, z, x) + T(z - 1, x, y) 
+ T(t(x - 1, y, z), t(y - 1, z, x), t(z - 1, x, y)) . 

The total number of expansions of the definition will then be 1 + 4T(x,y, z), because the latter 
function T(x, y, z) satisfies the recurrence 

T(x, y, z) = if x < y then 1 

else 1 +T(x - l,y,z) + T(y - l,z,x) + T(z - 1, x, y) 
+ T(t(x - 1, y, z), t(y - 1, z, x), t(z - 1, x, y)) . 

Before we analyze T(x, y, z) it will be helpful to consider a similar but simpler function 

V(x, y, z) = if x < y then 

else 1 + V(x - 1, y, z) + V(y - 1, z, x) + V(z - 1, x, y) . 

The function V(x, y, z) can be understood as follows. Construct a ternary tree by starting with 
a simple leaf containing the triple [x, y, z] and repeatedly applying the following operation: If any 
leaf [x, y, z] of the tree-so-far has x > y, attach the nodes 

[x~l,y,z], [y-l,z,x], [z-l,x,y] 

immediately below it. Then V(x, y, z) will be the number of nonleaf nodes in the final tree. (This 
function V(x, y, z) has been studied by Ilan Vardi [14]; some of his analysis is reproduced here.) 
[Note: I have an example in the MS., but I'll skip it unless I get goahead from Vlad later.] 
The evaluation of V(x,y, z) is trivial when x < y, and it's also fairly simple when x > y and 
x > z > y: In that case we have 

V(x,y,z) = l + V(x-l,y,z) = l + x-z + V(z-l,y,z). 

A further simplification arises when we realize that the values of V(x, y, z) are invariant if we 
translate all the parameters by any integer amount: 

V(x + l,y+l,z + l) = V(x,y,z). 
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Therefore we can shorten our notation and our discussion by assuming that min(x, y, z) = 0. 

Suppose the ternary tree has [x, y, z] at the root, where min(x, y, z) = and either x > y > z 
or z > x > y. Then all of its non-leaf nodes are of two kinds, 

A(a,b) = [a,b,0] 

„ . , , r , . where a > b > . 
B(a,b) = [6,0, a] 

Below A(a, b) are the three nodes 

[a- 1,6,0], [6-l,0,a], [-l,a,6] 

where [a — 1, 6, 0] is a leaf if a = 6 + 1, otherwise it is A(a — 1, 6); similarly [6 — 1, 0, a] is a leaf if 
6 = 1, otherwise it is B(a,b — 1); and [— I, a, 6] is always a leaf. Below B(a,b) are the same three 
nodes; they appear in a different order, but that does not matter. 

It follows that V(x,y, 0) = V(y, 0,x), for all x > y > 0, and that F(x,y,0) has a simple 
combinatorial interpretation: It is the number of lattice paths that start at {x, y) and stay within 
the set { (a, 6) | a > 6 > }. (A lattice path is a path in which each step decreases exactly one of 
the coordinates by unity.) We will say that such a lattice path is confined. 

Confined lattice paths can be enumerated by using Andre's well-known reflection principle 
(see, for example, [4, exerise 2.2.1-4]). Given x > y > y' > 0, the number of confined paths from 
(x,y) to a point (x',y') for some x' is equal to the number of all possible lattice paths from (x,y) 
to (y', y' — 1) minus the number of such paths that touch a diagonal point. Paths of the latter type 
are in one-to-one correspondence with lattice paths from (x,y) to (y' — l,y'); the correspondence 
is obtained by interchanging x moves with y moves after the diagonal is first encountered. Hence 
the number of confined paths from (x,y) to points of the form (x',y'), given x, y, and y', is 

x - y' + y - (y' - 1) \ ( x - {y' - 1) + y - y' 



x-y' J \ y-y 1 



and the total number of confined paths starting from (x, y) is 



E 

y'=l 



x + y + l-2y'\ ( x + y+l-2y' 



k=i 

This is the quantity V(x,y, 0). Notice that we have 



x-y' J \ y-y' 



x-y-l + 2k\ _ f x-y-l + 2k 
k { k-1 



vo. + ^-t((?)-(» 1 ))-i:(?) i TT-t*. 

k=i v v 7 v 77 k=i v 7 fc=i 

the sum of the first n Catalan numbers. 

Returning to the evaluation of T(x, y, z), we must add to V(x, y, z) the values T(t(a — 1, 6, 0), 
t(b—l, 0, a), a) at every node of type A(a, 6), and the values T(t(b— 1, 0, a), a, t(a — l, 6, 0)) at every 
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node of type B(a, 6). Fortunately these additional amounts are mostly zero. Our knowledge about 
the values of t(x, y, z) allows us to conclude that, when a > b > 0, we have 



t(a- 1,6,0) = a- 1; t(b -1,0, a) 



0, 6=1; 
a, b > 1. 



Therefore T(t(b — 1,0, a), a, t(a — 1,6,0)) = 0; and T(t(a— 1,6,0), t(b — 1,0, a), a) = except 
when 6 = 1. Only the nodes of type A(a, 1, 0) acquire additional values; and at such nodes we add 
T(a- 1,0, a). 

The number of non-root nodes of type A(a', 1) in a tree whose root is of type A(a, 6) or B(a, 6) 
is the number of nodes of type A(a' + 1, 1) or B(a' + 1, 1). And the number of such nodes is the 
number of confined lattice paths from (a, 6) to (a' + 1, 1), which is 

a - (a' + 1) + 6 - 1\ fa - 1 + 6 - (a' + 1) 
6-1 )~\ a-1 

by Andre's reflection principle. Therefore 

T(6, 0, a) = V(a, b, 0) + £ ( (° + ^ " ') " (° + " ') ) ^ ~ 1, 0, a') 



The same formula holds for T(a, 6,0), except that we must add T(a — 1,0, a) when 6 = 1; a root 
node of type A(a, 1) makes a contribution, but a root node of type B(a, 1) does not. Putting these 
facts together yields the following recurrence for the numbers T n = T(n,0,n + 1): 



T„ +1 ^,n + , n+1 ,o, + g ((»;*) -(:;:)) x„-, 



Let y ra = y(n + 1, n, 0); the first few values of these sequences are as follows: 

n=12345 6 7 8 9 
T4 = 1 3 8 22 64 196 625 2055 6917 
T n = 1 4 14 53 223 1034 5221 28437 165859 

It is not difficult to deduce that the numbers T n grow very rapidly, in fact faster than A n for any 
constant A: 

Lemma 2. If e > 0, we have 

T n > n^ n 

for all sufficiently large n. 

Proof. Choose k large enough so that k/(k + 1) > 1 — e. Looking only at the kth term of the 
recurrence for T n tells us that, for all n > k, we have 

ffn + k\ f n + k\\ _ (n + fc)(n + fe - 1) . (n + 2)(n + l- fe) 
T " +1> ii n J"ln + lJJ B_fc " ^ n " fe ' 



7 



Thus 

In T n+ i > ln(n + k) H h ln(n + 2) + ln(n + 1 - A;) - In A;! + In T„_ fc . 

Iterating this relation yields 

k 

lnT„ > nlnn + 0(n) 

k + 1 

and the result follows. (A similar but weaker result was obtained by Ilan Vardi [14], who used the 
fact that T n+ i > nT„_i to prove that lnT„ > ^nlnn + O(n).) □ 

In fact, we can prove a stronger lower bound by observing that for all n > 1, we have T n > b n , 
where b n is the Bell number defined by 



n—l ,\ 

b n+1 = 1 + ^ ( , J K-i, 
k=o ^ ' 



Each term in the recurrence for T n is greater than or equal to the corresponding term in the 
recurrence for b n , by induction. It is known [1, (6.2.7)] that 



i ^ n In n — n In In n — n 
Un ^ ^ 

for all sufficiently large n. Thus T n grows faster than (n/elnn) n . 
On the other hand, we can prove the upper bound 

T n < 3n! . 

This is clear for n < 3. Let t n = T n /nl, and rewrite the recurrence for T n as follows: 
V n +i 1 ( n + 2 

tn+l — 7 — 7T7 H -~r I t n + t n -\ H t n -2 



(n + 1)! n + 1 V 2n 

n + 3 n + 2 ra + 4 n + 3 n + 2 

+ 3(n - 1) ~2^~ tn ~ 3 + 4(n - 2) 3(n - 1) ~2^~ *"" 4 + 

The coefficients inside the parentheses are all < 1, so we have 

V n +i 1 3n Ki+i 

Wl < 7 1 \, + | 7 (tn + t n -l +-+tl)< | 7" + 



(n + 1)! n + 1 v " 17 n + 1 (n + 1)! 

by induction. And it is easy to verify that V n+ \ < 3n! for n > 3, because V n+ \ < 4™. (The Catalan 
numbers C n satisfy C n+ \ < 4C n , hence V n+ \ < C\ + 41/„ and we have V n+ \ < 4V n .) We have 
proved 

Theorem 2. When the Takeuchi recursion t(x, y, z) is used in a memoryless manner to evaluate 
t(n, 0, n + 1), the definition is expanded 1 + 4T(n, 0, n + 1) times, where 

e nlnn-n\nlnn-n < T ( n , Q, n + 1) < e ™l™-n+lnn 
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for all sufficiently large n. □ 

Only 0(n 2 ) evaluations are needed when previously computed results are remembered; thus memory 
is especially helpful here. 

A more precise asymptotic formula for V n can be obtained from the well-known generating 
function for Catalan numbers [3, page 203], 



C(z) 



n>l n>l 

1 - VI ~4:Z 



2z 

Darboux's lemma (see [5]) now shows that 



£ (H TV « • 

The generating function of the numbers T n satisfies a remarkable functional equation: We have 

T{z) = T n z n 

n 

= e v;, +l2 " +I = E (" t ") T ^ n+1 - E (I t f) T ^' +l 

n n,k ^ ' n,k ^ ' 

= + E (" V*) T » 2 " +H1 - E (; + _T) 



n,k n,k 

C{z) - 1 v n+1 C(^)" ^ ra+2 C(^+ 2 
l_ z + 2^* 2^"* ^T4I 



_ C(*)-1 ^( 2 -^)) TrcCL ^ 
because z(7(z) 2 = C{z) - 1 and £ fe ( n + 2fc ) z fc = C(z) V^4i. 

Open Problem 2. Obtain further information about the asymptotic properties of the coefficients 
T U T 2 ,.... □ 

The evaluation of t(x, y, z) turns out to be much, much faster if we apply the technique of 
lazy evaluation or "call by need" when expanding the definition. (See Vuillemin [15], [16].) Indeed, 
we can ignore the third argument t(z — l,x,y) in the recursion, unless we have discovered that 
t(x — l,y, z) > t(y — 1, z, x); so the number of times the else clause needs to be expanded satisfies 
the recursion 

K(x, y, z) = if x > y then 

else (l + K(x - 1, y, z) + K(y - 1, z, x) 

+ if t(x — l,y, z) < t(y — 1, z, x) then 

else K(z- l,x,y) +K(t(x- l,y,z), t(y - l,z,x), t(z - l,x,y))) . 
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And this recursion turns out to be quite simple. First, if x > y < z, we have 

K{x, y, z) = 1 + K(x - 1, y, z) = x - y 
because t(x — 1, y, z) < z and t(y — 1, z, x) = z. Second, if x > y > z + 1, we have 

K{x, y, z) = 1 + K(x - 1, y, z) + K(y - 1, z, x) 

= 1 + K(x - 1, y, z) + y - 1 - z = (x - y)(y - z) , 

because t(x — l,y,z) = x — 1 and t(y — 1, z, x) = x. Finally, if x > y = z + 1, we have i(x — 1, y, z) = 
x — 1, £(y — 1, z, x) = z, and t(z — 1, x, y) = x; hence 

K(x, y, z) = 1 + K(x — l,y, z) + x — 1 — z = (x — y)(x — y + 3)/2 . 

Incidentally, when x > y > z + 1, the expansions of else clauses occur at the arguments (£, y, z) 
and (77, z, £,) for x > £ > y and y > r) > z; when x > y = z + 1, they occur at (£, y, z) and (77, z, £) 
for x > £ > r/ > y. Since these arguments are distinct, no additional savings over call-by-need 
would be obtained by remembering previously computed values, unless t(x, y, z) is being evaluated 
at more than one point (x,y, z). The fact that the necessary arguments are limited underlies the 
simple mechanical proof of termination found by Moore [12]. 

3. False Takeuchi functions. Vardi [14] has considered a general recursion scheme of the form 

v h (x, y, z) = if x < y then h(x, y, z) else v h (v h (x - 1, y, z), v h (y - 1, z, x), v h (z - 1, x, y)) . 

If we set h(x,y,z) = 0, the function Vh(x,y,z) will of course be identically zero; we will deduce 
that its value is zero after expanding the definition exactly 1 + 4V(x, y, z) times, where V(x, y, z) is 
the function considered in the previous section. This is clearly the minimum number of expansions 
necessary over all possible auxiliary functions h(x,y,z), if we do not or cannot use call-by-need. 

Richard Gabriel used a Takeuchi- like function in extensive benchmark tests of Lisp compilers, 
but his function was slightly different from Takeuchi's original: 

g{x, y, z) = if x < y then z else g(g(x - 1, y, z), g(y - 1, z, x), g(z - 1, x, y)) . 

(Notice that in this case call-by-need is inapplicable.) Gabriel explains the discrepancy as follows 
[2, pp. 10-11]: 

When the Computer Science Department at Stanford University obtained the first two 
or three Xerox Dolphins, John McCarthy asked me to do a simple benchmark test with 
him. We sat down, and he tried to remember the Takeuchi function, which had had 
wide circulation. Because it was simple and because there were many results for it in the 
literature, he felt that it would be a good initial test. Of course, John misremembered the 
function. But we did not realize it until I had gathered a great many numbers for it. 

Indeed, Gabriel's book [2] gives detailed timings for the computation of y(18, 12, 6) = 7 on 132 differ- 
ent configurations, and he lists four additional variants of g that provide further types of benchmark 
tests. 
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The seemingly trivial change from t(x, y, z) to g(x, y, z) actually makes g(x, y, z) substantially 
easier to compute, if t(x, y, z) is not evaluated with memory of previous results or with call-by-need. 
Vardi [14] has shown that the corresponding running time G(n, 0, n + 1) is asymptotically less than 
(3 + y/8 ) n , although the exact order of growth is not known. Vardi has also observed that Gabriel's 
recursion defines the following curious pattern of values: 

g(x, y, z) = if x < y then z 

else if y > z then 

if y = z or (x — y) odd then y 
else z + 1 

else if z < x + 1 and (z < x or x > y + 1) then y 
else if (z — x) even then x 
else y + 1 . 

Here is another example of a generalized Takeuchi recurrence whose solution exhibits odd-even 
behavior: 

b(x, y,z) = if x < y then if x = y = z then else 1 

else b(b(x - 1, y, z), b(y - 1, z, x), b(z - 1, x, y)) . 

This time the output of the function is boolean — always either or 1 — although x, y, z range over 
all integers. The computed values turn out to be 

b(x, y,z) = if x < y then if x = y = z then else 1 

else if z > y + 1 then if (x — z) odd then else 1 
else if y = z then if (x — y) even then else 1 
else if (x — y) odd then else 1. 

The generalized recursion Vh does not always define a total function by repeated expansion. 
For example, consider the auxiliary function 

e(x, y, z) = if x odd then else 1; 

then we get 

Ue(l,0,0) =« e (v e (0,0,0), V e {-l,0,l), U e (-l,l,0)) =V e (l,0,0) 

and the recursion loops endlessly. There is a simple characterization of the cases where Vh is total 
in the boolean case: 

Lemma 4. Let h(x,y,z) map arbitrary integers x,y,z into or 1. Then the recursive equation 
for Vh(x, y, z) defines a total function Vh except in the following three cases: 

(i) h(0, 0, 0) = 1 and h(-l, 0, 1) = h(-l, 1, 0) = 0; 

(ii) h(0, 0, 1) = h(0, 1, 0) = 1 and h(-l, 1, 1) = 0; 

(iii) h(0, 0, 0) = h(0, 0, 1) = h(-l, 1, 0) = 1 and h(-l, 0, 1) = h(-l, 1, 1) = 0. 

Proof. If h(0, 0,0) = or h(— 1,0, 1) = 1, we have the well-defined result 

v h (l, 0, 0) = h(h(0, 0, 0), h(-l, 0, 1), h(-l, 1, 0)) ; 
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otherwise we have 

^(1,0,0) =^(1,0,^-1,1,0)) , 

which loops in case (i) but gives Vh(l, 0, 0) = ^(1, 0, 1) otherwise. Similarly, we find that h(0, 0, 1) = 
or h(— 1, 1, 1) = 1 implies 

v h (l, 0, 1) = h(h(0, 0, 1), h(-l, 1, 1, h(0, 1, 0)) ; 

hence ^(1,0, 1) is well defined whenever case (ii) does not hold, except in the case when it leads 
to ^(1,0,1) =^(1,0,0). 

If neither case (i) nor case (ii) holds, then Vh{x,y,z) is well defined for all boolean values 
x,y,z, except in case (iii). And when Vh(x,y,z) is defined for all boolean x,y,z, we can evaluate 
Vh(x, y, z) for all x, y, z in 0(V(x, y, z)) steps. □ 

When the boolean function Vh of Lemma 3 isn't total, we can always complete it to a total 
function Vh(x,y,z) that does satisfy the recurrence. We simply assign arbitrary boolean values to 
^(1,0,0) and/or ^(1,0,1), whichever is undefined. For example, there are four total functions 
v e (x, y, z) that satisfy the recurrence arising from the auxiliary function e(x, y, z) considered above: 

(x, y, z) = if x odd then else 1; 

v^ 1 (x,y,z) = if x < y then if x odd then else 1 
else if x even then 1 
else if y odd then 
else if z odd then 1 else 0; 

i>g°(x, y, z) = if x < y then if x odd then else 1 
else if x odd then 

if y odd or z odd then else 1 
else if y odd or z odd then 1 
else if y < z < x then 1 else 0; 

vl 1 (x,y,z) = if x < y then if x odd then else 1 
else if x even then 1 
else if y odd then 
else if z < y then 1 
else if z odd then 1 else 0. 
However, the recurrence Vh(x,y,z) cannot be completed to a total function for arbitrary aux- 
iliary functions h(x,y,z). Consider, for example, the (admittedly contrived) mapping 

h(x, y, z) = 2xy - 4x + y + z — 1 . 

There is no total function Vh because we would otherwise have 

« h (2,l,4) =v h (h(l, 1,4), /i(0,4,2), v h (h(2,2,l), ^(1,1,3), h(0,3,2))) 
= ^(2,5,^(2,1,4)) = 16 + ^(2,1,4). 



12 



An incompletable function can even be constructed when we restrict ourselves to auxiliary 
functions that are limited by the condition 

h(x, y, z) < max(x, y, z) . 

For example, suppose we define 

h(x, y, z) = if (x, y, z) = (1, 1, 4) then 4 

else if (x, y, z) = (3, 3, 3) then 2 
else if (x, y) = (2, 3) then 1 
else if max(x, y, z) > 3 then 3 
else max(x, y, z) . 

Then we have Vh(x, y, 3) =3 whenever 3 > x > y. For if x = y, clearly Vh(y, y, 3) = h(y, y, 3) = 3; 
otherwise 

v h (x,y,3) = v h (v h (x - l,y,3), v h (y - 1,3, x), v h (2,x,y)) 
= v h (3,3,< 2) = 3. 

Therefore if y < 3 we have 

^(3,2/, 3) = v h (v h (2,y,3), v h (y - 1,3,3), ^(2,3,y)) 

= ^(3,3, 1) = 3. 

It follows that 

v h (3,y,3) = if y = 3 then 2 else 3 . 

But we also must have 

^(4,3,1) =v h (v h (3,3,l), ^(2,1,4), ^(0,4,3)) 

= v h {3, v h (v h (l,l,4), ^(0,4,2), ^(3, 2,1)), 3) 

= v h (3, ^(4,3,^(^(2,2,1), ^(1,1,3), ^(0,3,2))), 3) 

= v h (3, ^(4,3,^(2,3,3)), 3) 

= v h (3, ^(4,3,1), 3) . 
And there is no y such that y = Vh(3, y, 3). 

Open Problem 3. If we restrict h(x,y,z) to be strictly less than max(x,y,z), is there always a 
total function Vh(x,y,z) that satisfies the generalized Takeuchi recurrence? □ 

We have considered Takeuchi's special case h(x, y,z) = y as well as Gabriel's special case 
h(x, y, z) = z, so it is natural to consider also the recurrence with h(x, y, z) = x. Let 

k(x, y, z) = if x < y then x 

else k(k{x — 1, y, z), k(y — 1, z, x), k(z — 1, x, y)) . 

This recursive definition yields only a partial function because, for example, we have 

k(x + 1, x, x) = k(x, x — 1, x — 1) = k{x — 1, x — 2, x — 2) = • • • . 
However, there are infinitely many ways to define a total function that does satisfy the k recurrence: 
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Theorem 3. Let c be any integer. The function 

k c (x, y, z) = if x < y then x 

else if y < z + 1 then c else min(y, c) 

satisfies the generalized Takeuchi recurrence stated above for k(x,y, z). 
Proof. Notice that we have the special values 

k c (x, c, z) = min(x, c) ; 

if x > y then k c (x, y,c) = c . 

The proof is now by induction on x — y. 
If x = y + 1 we have 

k c (k c (x - l,y,z), k c (y - l,z,x), k c (z - l,x,y)) = k c (y, k c (y -l,z,y + 1), k c (z -l,y+ l,y)) . 
If y < z + 1, this reduces to 

k c {y,y ~ l,z- 1 or c) = c; 

and if y > z + 2, it is 

k c (y, c, k c (z -l,y+ l,y)) = min(y,c) . 

Thus we obtain k c (x, y, z) when x = y + 1. 
If x > y + 2 and y < z + 1 we have 

k c (k c (x- l,y,z), k c (y - l,z,x), k c (z - l,x,y)) = k c (c, y-1, z- 1 or cor x) , 

which equals c since y — \ < z — 1 + 1 and y — 1 < x + 1 . 

And finally if x > y + 2 and y > z + 2 the right side of the recurrence reduces to 

k c (min(y, c),c,z-l) = mm(y, c) . □ 

Corollary. The least fixed point of the recursive definition k(x, y, z) is if x < y then x else lo. 
Proof. Whenever x > y, we have k c (x, y,z) = y when c = y but not when c < y. □ 

4. The Takeuchi recurrence in higher dimensions. If we define 

t(w, x, y, z) = if w < x then x 

else t(t(w — 1, x, y, z), t(x — 1, y, z, t(y — 1, z, w, x), t(z — 1, x, y)) 

it turns out that the function reduces to the simple mapping 

t(w, x, y, z) = if w < x then x else if x < y then y 
else if y < z then z else u; . 
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Therefore it is natural to conjecture that the m-dimensional generalization 

t(x\, x 2 , ■ ■ ■ , x m ) = if x\ < x 2 then x 2 

else t(t(x 1 - l,x 2 , ■ ■ -,x m ), t(x m - l,xi, . . . ,x m _i)) 

is satisfied by the m-dimensional "first rise" function 

u(xi,X2, • • • , x m ) = if X\ > ■ ■ ■ > x k < Xk+i for some k > 1 then Xk+i 
else x\ . 

But this is false, for all m > 4. Indeed, we have a 5-dimensional counterexample, 
t(5, 3, 2, 0, 1) = i(t(4, 3, 2, 0, 1), 2, t(l, 0, 1, 5, 3), 1, 5) 

= t(t(3,2,t(l,0,l,4,3),l,4),2, t(0,l,...), 1, 5) 

= i(t(3,2,t(0,l,...),M), 2, 1, 1, 5) 

= t(t(3,2,l,l,4), 2, 1, 1, 5) 

= t(t(2, 1,1,4,3), 2, 1, 1, 5) 

= t(t(l, 1,4,3,2), 2, 1, 1, 5) = 4(1,2,1,1,5) =2, 

while u(5,3,2,0,l) = 1. 

The true general behavior is somewhat complicated, although (fortunately) the complications 
do not get worse and worse as m grows larger and larger. Let us define an auxiliary set of functions 
gj(x±, . . . , Xj) for j > 2 as follows: 

gj (x\ , . . . , Xj) = if j = 2 then x 2 

else if xi = x 2 + 1 then gj_\{x2, ■ ■ ■ ,Xj) 
else if x 2 = x 3 + 1 then max(x 3 ,Xj) else Xj . 

Theorem 4. The function 

f(xi, . . . , x m ) = if xi > ■ ■ ■ > Xk < for some A; > 1 
then g k+1 (xi, . . . ,x k +i) else x\ 

satisfies then m-dimensional Takeuchi recurrence. 

Proof. Given x\, . . . , x m , with x\ > x 2 , let 

Dj fip^j ^ ' ^J + l > • • • i •^m j 2-1 ; • • • j •Cj — 1 ) • 

We want to show that f(y 1 , . . . ,y m ) = f(x 1 ,. . .,x m ). 

If Xi > ■ ■ ■ > x m , we have yi = x\ — 1 or x 2 ; y 2 = x\ or x 3 ; . . . ; y m _i = ari or x m ; and 
Vm = x\. We cannot have yi > • • • > y m , because y m > yi- Hence there is a unique k > 1 such 
that yi > ■ ■ ■ > yfc < j/fe+i- And this can happen only if y^+i = ^l! otherwise y k +i = x k+2 < y k . 
It follows that /(yi, . .. ,y m ) = g k +i(yi, ■ • ■ ,Vk+i) = x 1 = f(x u ...,x m ). 
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Assume now that x\ > • • • > x k < x k +i, where k > 2, and let a be as large as possible such 
that Xi = Xi+\ + 1 for 1 < i < a. Then y, = Xi+i for 1 < i < a. If a = k, we have y k = x k +i and 
f(yi, ■ ■ -iVm) = 9k(yi,- ■ ■ ,2/fc) = 9k(x 2 , ■ ■ -,Xk+i) = g k +i(xi,. . .,x k +i) = /(xi, • • • ,x m ). 

Assume therefore that a < k; hence x a > x a +i + 1. Let b be as large as possible such that 
Xi > Xi+i + 1 for a < i < b. If b = k, we have yi = Xk+i for a < i < b, hence f(yi, ■ ■ ■ ,y m ) = 

9a+l(yi, ■ ■ ■ , 2/o+l) = 9a+l(x 2 , . . . ,X a ,X k+1 ,X k+1 ) = X k+i = g k+1 (x 1 , . . . ,Xfc+l) = f(xi, . . .,x m ). 

Assume therefore that b < k; hence Xb = + 1 and y^ = x^+i- Let z = max(xf,+i, x k+ \). 
We have = z for a < i < b - 1. If a > 1 and z > x a , we have f(yi, ... , y m ) = g a (vi, ■ ■ ■ ,y a ) = 
g a {x 2 ,.. .,x a ,z)=z = g k +i(xi,. . -,x k+1 ) = f{x u . . .,x m ). 

Assume therefore that a = 1 or z < x a . If Xb-i > Xb + 2 then yb-\ = z; otherwise yb-i = 
xi > x k+ i for some I in the range 6+l<Z<fc + l. If 6 > a + 2, we have f(y\, ■ ■ ■ ,y m ) = 
g a +i(yi, ■ ■ -.Z/o+i) = 9a+i(x2, ... ,x a ,z,z) = z = g k+ i{xi, . . .,x k+1 ) = f(xt, . . . ,x m ). If b = a + 2, 
the same chain of equalities is valid unless Xb-\ = Xb + 2 and yb-i = x\ and xi 7^ z. In the latter 
case we cannot have z = x k+ i, for that would imply x k+ \ > Xb+i > xi, hence xi = x k +i = z. 
It follows that z = x b+ i > x k+1 , and x x < z = y b . Then f(yi,...,y m ) = 50+2(2/1, 2/0+2) = 
g a+2 {x 2 ,.. .,x a ,z,x h z) = z = f(xi, . . .,x m ). 

Assume therefore that b = a + 1. If z = x k +i we have x k+ \ > Xi for b < i < k, hence y a = z 
and yi = Xi + \ or z for a < i < k. The first appearance of z among y a +i, • • • , 2/fc+i will show that 
/(2/1, ...,y m ) = z = /(xi, . . . ,x m ). 

Assume therefore that z = Xb+i > x k+ \. If x a > Xb+2 then y a = yb = z, hence f(yi, . ■ ■ , y m ) = 
g b (x 2 , ... ,x a ,z,z) = z = /(xi, . . . ,x m ). 

Assume therefore (and finally) that x a = x^ + 2, so that y a = X\ > x k+ i, where 6+1 < I < k+1. 
Then xi < x b+ i = yb, and /(yi, . . . ,y m ) = g b (x 2 , . . . , x a , x u x b+1 ) = z = /(xi, . . . ,x m ). We have 
proved that f(yi, ■ ■ ■ , y m ) = f(xi, ■ ■ ■ , x m ) in all cases. □ 

A machine-based proof of Theorem 4 would be very interesting, especially if it could cope with 
functions having a variable number of arguments. 

Notice that we have not proved that the m-dimensional Takeuchi recursion t(x\, . . . ,x m ) ac- 
tually defines a total function, when m > 3. We have only shown that /(xi, . . . ,x m ) satisfies 
the recurrence. If the repeated expansion of t(x\, . . . ,x m ) actually terminates for some sequence 
of arguments (x±, . . . ,x m ), it must yield the value f(xi, . . . ,x m )', but we have not demonstrated 
that termination will occur, and there is apparently no obvious ordering on the integer m-tuples 
(xi, . . . , x m ) that will yield such a proof. Therefore we come to a final question, which will perhaps 
prove to be the most interesting aspect of the present investigation. 

Open Problem 4. Does the m-dimensional Takeuchi recursion equation define a total function, 
for all m > 3, if it is expanded fully (without call- by- need)? Equivalently, does the recurrence 
T(x\, . . . , x m ) = if Xi < x 2 then 

else 1 + T(x 1 - l,x 2 , . . . ,x m ) +T(x 2 - 1, x 3 , . . . , x m , xi) H 

+ T(x m —i 1, x m , Xi, . . . , x m — 2 ) + T[x m 1, X\, . . . , x m —i) 
+ T(f{x 1 - l,x 2 ,...,x m ),...,/(x m - l,xi,...,x m _i)) 
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define a total function on the integers (xi, . . . , x m ), for all m > 3? (Here / is the function of 
Theorem 4.) □ 

Close inspection of the proof of Theorem 4 implies that a call-by-need technique will always 
terminate when applied to the recursive equation for t(xi, . . . , x m ). If x\ > x 2 > • • • > Xk < Xk+i, 
the values yi = t(xi — . . . need be expanded only for 1 < i < k + 1, and this will 

be sufficient to determine the value of t(y 1 , . . . , y m ) = t(x\, . . . ,x m ) in a finite number of steps. 
(The proof is by induction on k.) However, the possibility remains that an attempt to expand the 
"irrelevant" parameters yk+2, ■ ■ ■ ,Um might loop forever. If so, the Takeuchi recurrence would be 
an extremely interesting example to include in all textbooks about recursion. 
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